Dark Matter in the Singlet Extension of MSSM: 
Explanation of Pamela and Implication on Higgs Phenomenology 

Wenyu Wang 1 , Zhaohua Xiong 1 , Jin Min Yang 2 , Li-Xin Yu 2 
1 Institute of Theoretical Physics, College of Applied Science, 
Beijing University of Technology, Beijing 100020, China 
2 Key Laboratory of Frontiers in Theoretical Physics, 
Institute of Theoretical Physics, Academia Sinica, Beijing 100190, China 



Abstract 

As discussed recently by Hooper and Tait, the singlino-like dark matter in the Minimal Su- 
persymmetric Standard Model (MSSM) extended by a singlet Higgs superfield can give a perfect 
explanation for both the relic density and the Pamela result through the Sommerfeld-enhanced 
annihilation into singlet Higgs bosons (a or h followed by h — > aa) with a being light enough to 
decay dominantly to muons or electrons. In this work we analyze the parameter space required 
by such a dark matter explanation and also consider the constraints from the LEP experiments. 
We find that although the light singlet Higgs bosons have small mixings with the Higgs doublets 
in the allowed parameter space, their couplings with the SM-like Higgs boson hsM (the lightest 
doublet-dominant Higgs boson) can be enhanced by the soft parameter A K and, in order to meet 
the stringent LEP constraints, the hsM tends to decay into the singlet Higgs pairs aa or hh in- 
stead of bb. So the hsM produced at the LHC will give a multi-muon signal, h$M —* aa — > 4/x or 
hsM — > hh — > 4a — > 8fi. 

PACS numbers: 14.80.Ly,11.30.Pb,95.35.+d 



I. INTRODUCTION 



The experiment Pamela has observed an excess of the cosmic ray positron in the energy 
range 10-100 GeV [1], which is hard to be explained by the conventional cosmic ray source 
[2]. While there may exist some mundane explanations like pulsars [3] and the acceleration 
of positron secondaries in cosmic ray acceleration regions [4] , the dark matter interpretation 
[5, 6] is especially interesting since it may be related to new physics to be probed at the 
LHC. 

To explain the Pamela excess by the dark matter annihilations, there are some challenges. 
First, the dark matter must annihilate dominantly into leptons since Pamela has observed 
no excess of anti-protons [1] (However, as pointed in [7], this statement may be not so solid 
due to the significant astrophysical uncertainties associated with their propagation). Second, 
the explanation of Pamela excess requires an annihilation rate which is too large to explain 
the relic abundance if the dark matter is produced thermally in the early universe. To tackle 
these difficulties, a new theory of dark matter was proposed in [6]. In this new theory the 
Sommerfeld effect of a new force in the dark sector can greatly enhance the annihilation 
rate when the velocity of dark matter is much smaller than the velocity at freeze-out in the 
early universe, and the dark matter annihilates into light particles which are kinematically 
allowed to decay to muons or electrons. 

The above fancy idea is hard to realize in the popular Minimal Supersymmetric Standard 
Model (MSSM) because there is not a new force in the neutralino dark matter sector to 
induce the Sommerfeld enhancement and the neutralino drak matter annihilates largely to 
final states consisting of heavy quarks or gauge and/or Higgs bosons [8, 9]. However, as 
discussed in [10], in the extension of the MSSM by introducing a singlet Higgs superfield, 
the idea in [6] can be realized by the singlino-like neutralino dark matter (hereafter the 
singlino-like neutralino is simply called singlino): 

(i) The singlino dark matter annihilates to the light singlet Higgs bosons and the relic 
density can be naturally obtained from the interaction between singlino and singlet 
Higgs bosons; 

(ii) The singlet Higgs bosons, not related to electroweak symmetry breaking, can be light 
enough to be kinematically allowed to decay dominantly into muons or electrons 
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through the tiny mixings with the Higgs doublets; 

(iii) The Sommerfeld enhancement needed in the dark matter annihilation for the expla- 
nation of Pamela result can be induced by the light singlet Higgs boson (h). 

Such an explanation of dark matter requires that the singlet Higgs field has very small 
mixing with the Higgs doublets, which implies that the singlino dark matter may remain 
hidden and irrelevant to the LHC experiments. However, we note that the singlet extension 
of the MSSM has a quite large parameter space and thus the coupling of the light singlet 
Higgs (h, a) with the doublet Higgs (the lightest one is called hsn) may be enhanced by 
other parameters. For example, through the soft term A K S 3 (S is the singlet Higgs field) 
with a large A K , a pair of singlet Higgs bosons may sizably couple to a doublet Higgs boson 
although the mixing between the singlet and doublet Higgs fields is small. Therefore, this 
model may allow for exotic Higgs phenomenology at the LHC. 

In this work we study the parameter space allowed by the explanation of Pamela result 
plus relic density via Sommerfeld enhancement and also consider the constraints from the 
LEP experiments. We find that although the light singlet Higgs bosons have small mixings 
with the Higgs doublets, their couplings with the SM-like Higgs boson (hsu) can be enhanced 
by the soft parameter A K and, in order to meet the stringent LEP constraints, the hsM 
tends to decay into the singlet Higgs pairs aa or hh instead of bb. This implies that the hsM 
produced at the LHC will give a multi-muon signal, hsM da — > 4/x or Hsm —^hh-^Aa^ 
8/i. 

This work is organized as follows. In Sec. II we discuss the Higgs and neutralino sectors 
in the singlet extension of the MSSM. In Sec. Ill we scan the parameter space allowed by 
the dark matter explanation and LEP experiments, and discuss the implication on Higgs 
phenomenology. Finally, a summery is given in Sec. IV. 

II. HIGGS AND NEUTRALINOS IN SINGLET EXTENTION OF MSSM 

The Higgs superpotential in the general singlet extension of the MSSM is given by [10] 

W = i^H u -H d + XSH u -H d + r ] S+^fi s S 2 + ^KS 3 , (1) 

where S is the singlet Higgs superfield while H u and H d are the doublet Higgs superfields. 
The Higgs scalar potential consists of the D-term, the F-term and the soft SUSY-breaking 
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term. Since S is a singlet, the D-term is same as in the MSSM. The F-term from the 
superpotential is given by 

V F =\fi + XS\ 2 (\H U \ 2 + \H d \ 2 ) + + fi s S + XH U ■ H d + kS 2 \ 2 . (2) 

The soft SUSY-breaking terms are given by 

V so {t = m 2 Ru \H u \ 2 + ml d \H d \ 2 + m 2 s \S\ 2 

+ (BfiH u ■ H d + \A\ H u ■ H d S + CrjS + X - B sf i s S 2 + ^ kA k S 3 + h.c.) . (3) 

So the Higgs potential reads 

V = \ii + \S\ 2 (\H u \ 2 + \H d \ 2 ) + \\H u -H d + KS 2 \ 2 

+m 2 iu \H u \ 2 + m 2 ld \H d \ 2 + ml\S\ 2 

+ (BfiH u ■ H d + XA X H U ■ H d S + Xfi s H u ■ H d S* 

+C V S + - B sf i s S 2 + - kA r S 3 + Kfi s S 2 S* + h.c.) (4) 

where g 2 = (gf + g%)/2 with g 1 and g 2 being respectively the coupling constant of SU(2) and 
U(l) in the SM. 

After the Higgs fields develop the vevs h u , h d and s, i.e., 

rro , . H uR + iH uI H dR + iH dI S R + iSi 

K = h u + -j= , H» = h d + — , S = s + — -j=— (5) 

we obtain a 3 x 3 mass matrix Aih for CP-even Higgs bosons, a 3 x 3 mass matrix Ai a for 
CP-odd Higgs bosons and a 2 x 2 mass matrix M. c for the charged Higgs bosons: 

(1) The CP-even Higgs mass matrix in the basis (H uR , H dR , S R ) is given by 

M h , u = g 2 h 2 u + cot (3 [Xs(A x + ks + fi s ) + B/j] , (6) 

M h)2 2 = g 2 h 2 d + tan/5 [Xs(A x + ks + fi s ) + Bfx] , (7) 

M h , 33 = KA X + /is )^ - \»(h 2 u + h 2 d ) + ks{A k + Aks + 3/i s ) - ^5, (8) 

s s s 

M K12 = (2A 2 - g 2 )h u h d - Xs{A x + ks + fi s ) - Bp, (9) 
M h , 13 = 2X(/2 + Xs)h u - Xh d (A x + 2ks + fi s ), (10) 
M h , 23 = 2A(/i + Xs)h d - Xh u (A x + 2ks + fi s ), (11) 



where tan/3 = h u jh d . This mass matrix can be diagonalized by a rotation 



U 



HdR 



(12) 



with an orthogonal matrix U. The mass eigenstates are ordered as < < m^ 3 . 
In the MSSM limit (A, rj, /i s , k — >• and /i 3 ~ S^) the elements of the first 2x2 
sub-matrix of U are related to the MSSM angle a as 



cos a , 
— sin a , 



U 2 i = sin a , 
f/ 2 2 = cos a . 



(2) The CP-odd Higgs mass matrix M. a in the basis (H u i, Hji, Si) is given by 
M a ,ii = cot (3[\s(A x + KS + fi s ) +B(l], 

M a ,22 = t&nj3[\s(A x + KS + Us) + Bfj], 

h h 

M a , 33 = A\Kh u h d + \(A X + /i s )^^ 

s 

-X-(h 2 u + h 2 d ) - ks(3A k +n 8 )-?!- 2B s fi s , 
s s 

M a ,12 = \s(A X + K,S + fi s ) + Bfl, 

M a ,i3 = Xh d (A x -2Ks - n s ), 

M a ,23 = \h u (A X -2KS - ll s ). 



(13) 

(14) 
(15) 

(16) 
(17) 
(18) 
(19) 



The diagonalization of this mass matrix can be performed in two steps. The first step 
is to rotates into a basis (A, G, Si) with G being a massless Goldstone mode: 



cos/3 —sin/? 
sin p cos (5 
1 



G 
\ Sl J 



(20) 



H u i 
Hdi 

\ Sl ) 

Dropping the Goldstone mode, the remaining 2x2 mass matrix in the basis (A, Si) 
is given by 



MnU 



(tan P + cot P) [\s(A x + ks + fi s ) + Bfj] , 



M a ,22 = 4:\Kh u h d + X(A X + n s ) 

C V 



h u h 



-ks(3A k + /jl s ) 2B s fi a , 



M a ,i2 = Xyjhl + h? d (A x -2ks- n s ). 



(21) 

(22) 
(23) 



It can be diagonalized by an orthogonal 2x2 matrix P' and the physical CP-odd 

states cij are given by (ordered as m ai < m a , 2 ) 

ai = P[ X ~A + P[ 2 Sj = /^(cos/WZuj + sin (3H dI ) + P[ 2 Sj, (24) 
a 2 = P' 21 A + P' 22 S! = P' 21 (cos !3H uI + sin f3H dI ) + P^Sj, (25) 

(3) The charged Higgs mass matrix M± in the basis (if if j~) is given by 

M±= (\s(A x + Ks + fi s )+Bv + h u h d $-\ 2 )) \ COtP ' ] (26) 
V 2 / \ 1 tan/3 



which gives one eigenstate H ± of mass Tr.M± and one massless goldstone mode G 



if* = cos - sin f3G ± , 



ff± = sin (3H± + 008/36^ . 



(27) 



(4) The neutralino mass matrix A4o can be read from the Lagrangian 



1 loo 

- M1A1A1 + - M2A2A2 



(28) 



where Ai is the U(l)y gaugino and A2 is the neutral £77(2) gaugino. In the basis 
ip° = (— — iX 2 , ipui V's) we obtains 

1 



£ = -^Wf + h.c, 



(29) 



where 



-Mo 



Mo 



\ 



gih u 

V2 


gihd 


\ 


92 h u 


92h d 
V2 










-Xh d 







-Xh u 
2ks + jj, s J 



(30) 



Diagonalizing this mass matrix, one obtains 5 mass eigenstates (ordered in mass) 

X° = (31) 
6 



III. EXPLANATION OF PAMELA AND IMPLICATION ON HIGGS DECAYS 



In our study the lightest CP-odd neutral Higgs boson a\ is singlet-dominant, while for the 
CP-even neutral Higgs bosons the lightest one hi is singlet-dominant and the next-to-lightest 
h,2 is doublet-dominant. We use the notation: 



eti, 



h = hi, 



h 



SM 



(32) 



As discussed in [10], when the lightest neutralino Xi i n Eq.(31) is singlino- dominant, it can 
be a perfect candidate for the dark matter. As shown in Fig.l, such singlino dark matter 
annihilates to a pair of light singlet Higgs bosons followed by the decay h — > aa (h has very 
small mixing with the Higgs doublets and thus has very small couplings to the fermions). In 
order to decay dominantly into muons, a must be light enough. Further, in order to induce 
the Sommerfeld enhancement, h must also be light enough. From the superpotential term 
kS 3 we know that the couplings hxiXi an d aXiXi are proportional to k. To obtain the relic 
density of the dark matter, k should be 0(1). 



X? FT* 



h, a 
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FIG. 1: Feynman diagrams for singlino dark matter annihilation where Sommerfeld enhancement 
is induced by exchanging h. 

Since h, a must be singlet-dominant and Xi must be singlino- dominant, this implies small 
mixing between singlet and doublet Higgs fields. From the superpotential in Eq.(l) we see 
that this means the mixing parameter A must be small enough. On the other hand, the 
smallness of A is also required by the lightness of hi and ai whose masses are approximately 
given by 



M 



h,33 



M 



a,22 



KS 



KS 



A(^A + fi 



h,,,h 



u lb d 



^-^(h't + h 2 d )+(A K + Aks + 3/i s - ^ 



A(^A + fi 



s 2 
KS Z 

h u hd 

s) T~ 

KS A 



fi 



~ ^(K + K) - [3A K + fi s + —L + 



(33) 
(34) 



KS* \ KS* KS 

In the following we scan over the parameter space. We modify the package NMSSMTools 
[11] and use it in our calculations. As discussed above, A must be small enough in order to 



get a singlino-dominant Xi an d singlet-dominant h, a (we checked from our scan that A must 
be smaller than 0.01 in order to get m a < 0.5 GeV and m h < 20 GeV). So in our following 
scan we fix A = 10~ 3 . Further, k is taken as 0.5, and for the squark sector the soft masses 
and the trilinear terms are fixed as 500 GeV. Other parameters vary in the ranges: 

-500 GeV < C, fi, Us, B, A x , M u M 2 < 500 GeV 

-(500 GeV) 2 < r] < (500 GeV) 2 , s < 500 GeV, 2 < tan/3 < 40. (35) 

In order to get small M.n,zz and AA a ,22, the third terms in Eqs. (33,34), which are not sup- 
pressed by a small A, must also be small. Therefore, in our scan we require parameters A K 
and B s to be in the ranges: 

A K e (-Aks - 3u s + ^X) ± 20GeV, (36) 

\ KS Z J 

2B s fi s E (^-3A k ks - fi s Ks - ± (3GeV) 2 (37) 
In addition, we consider the following constraints: 

(i) The constraints from the LEP experiments, which include the LEP1 bound on invisible 
Z decay and the LEP2 direct searches for Higgs bosons; 

(ii) m ai < 0.5 GeV; 

(iii) The singlino-like Xi to give the dark matter relic density ^oh 2 in the range 0.01-0.2, 
which can be calculated from the approximate formula [10] 

^--(^(aW) 2 - <w 

To calculate the Sommerfeld enhancement we follow [6] to numerically solve the Schrodinger 
equation 

! ' X + V(r) X =—X (39) 



2.1/ dr- x v 2M 
with the boundary condition (r — > 00) 

x(r) — > sin(/cr + 5), (40) 

where M and k are respectively the mass and momentum of the dark matter particle. V(r) 
is the Yukawa potential induced by exchanging h and is given by 



The Sommerfeld enhancement is then given by 



T 



(42) 



The survived points are displayed in different planes in Figs. 2-6. We see from Fig. 2 that 
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FIG. 2: The scatter plots showing the decay branching ratios a — * (muon), a — > gg (gluon) 

and a — > ss (s-quark) versus m a for A = 10 -3 . 



in the range 2m ^ < m a < 2m^-, a decays dominantly into muons. From Fig. 3 it is clear that 
h can be as light as a few GeV, which is light enough to induce the necessary Sommerfeld 
enhancement as shown in Fig. 4. In the calculation of the Sommerfeld enhancement, we 
assumed the dark matter move with a velocity 150 km/s. 

The fit to Pamela result has been given in [10]. As shown in Table I in [10], for the 
parameter space in Figs. 2-4 with 2m M < m a < 2m n and as light as a few GeV (so 
the Sommerfeld enhancement factor is large enough), the Pamela positron excess can be 
naturally explained. 

In Fig. 5 we show the branching ratios of hsM decays. We see that in the allowed parameter 
space hsM tends to decay into aa or hh instead of bb. This can be understood as following. 
The MSSM parameter space is stringently constrained by the LEP experiments if Hsm is 
relatively light and decays dominantly to bb, and to escape such stringent constraints hsM 
tends to have exotic decays into aa or hh. As a result, the allowed parameter space tends to 
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FIG. 3: Same as Fig.2, but showing and rrih„ versus m a 
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FIG. 4: Same as Fig.2, but showing the Sommerfeld enhancement factor induced by h. 



favor a large A K , as shown in Fig.6, which greatly enhances the couplings hsM^a and hsivihh 
through the soft term kA k S 3 although S has a small mixing with the doublet Higgs bosons. 
Such an enhancement can be easily seen. Take the coupling hsMhh as an example. The soft 
term kA k S 3 gives a term kA k S r which then gives the interaction kA k L^L^ hsMhh because 
Sr = Uishi + + ^33^3 with hi = h and h 2 = hsM (see Eqs.12 and 32). Although the 



10 




100 102 104 106 108 110 112 



m^JGeV) 

FIG. 5: Same as Fig.2, but showing the branching ratio of h$M decays. The V (blue), 'x' (green) 
and '+' (red) denote the branching ratios of hgM — > aa, hsM — ► hh and hsM — ► bb, respectively. 




0.5 0.6 0.7 0.8 0.9 



Br(h SM -^aa, hh) 

FIG. 6: Same as Fig.2, but showing \A K \ versus the branching ratio of hsM aa,hh. 

mixing Uf 3 U23 is small for a small A, a large A K can enhance the coupling hshdhh. 

The SM-like Higgs boson hsM will be intensively searched at the LHC and its dominant 
decay mode in the MSSM is bb. In the singlet extension of the MSSM, its dominant decay 
mode may be changed to aa or hh, as shown in our above results. Such new decay modes 
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will give a multi-muon signal for Hsm at the LHC, i.e., Hsm — > clcl — > 4/i or hsM hh ^ 
Aa — > 8/1. So the phenomenology of h S u will be quite different from the MSSM predictions. 
Finally, we make some remarks regarding our results: 

(1) The recent DO search for h — > aa — > 4/i or 2/i2r channel obtained null results, which 
constrained the parameter space for the CP-odd Higgs a in the mass range of 3.6-9.5 
GeV [12]. But they do not constrain the parameter space considered in our analysis 
because we considered a much lighter CP-odd Higgs a with a mass below 0.5 GeV. 
Also, as pointed in [10], such a light a is allowed by T(3s) — > 7a — > / ~ffi + fi~ [13] and 
fT + — > 7r + a — > 7r + /x + /x~ [14] because in our scenario a is over dominated by singlet. 

(2) In the allowed parameter space displayed in our results, the mass of the SM-like Higgs 
boson hsM is rather below its theoretical upper bound (about 135 GeV in the MSSM). 
The reason is that, in order to push up its mass, the loop effects of heavy stops are 
needed ( note that in the singlet extension the tree-level upper bound can be enhanced 
by a term proportional to A, which is very small in our scenario). In our calculations 
the soft mass parameters in the squark sector are fixed to be 500 GeV and hence the 
stops are not heavy enough to push the mass of Hsm up to 135 GeV. Of course, we can 
choose heavy stops to push up the mass of hsM, in which case the allowed parameter 
space displayed in our results (with a relatively light h SM decaying dominantly into 
aa or hh) can still survive. 

(3) For the specified singlet extensions like nMSSM and NMSSM [15], the explanation 
of Pamela and relic density through Sommerfeld enhancement is not possible. The 
reason is that the parameter space of such models is stringently constrained by various 
experiments and dark matter relic density [16], and, as a result, the neutralino dark 
matter may explain either the relic density or Pamela, but impossible to explain both 
via Sommerfeld enhancement [17]. For example, in the nMSSM various experiments 
and dark matter relic density constrain the neutralino dark matter particle in a narrow 
mass range [16], which is too light to explain Pamela. 
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IV. SUMMARY 



The singlino-like dark matter in the MSSM extended by a singlet Higgs superfield can 
give a perfect explanation for both the relic density and the Pamela result through the 
Sommerfeld-enhanced annihilation into singlet Higgs bosons (a or h followed by h — > aa) with 
a being light enough to decay dominantly to muons. In this work we analyzed the parameter 
space allowed by such a dark matter explanation and also considered the constraints from 
the LEP experiments. We found that although the light singlet Higgs bosons have small 
mixings with the Higgs doublets in the allowed parameter space, their couplings with the 
SM-like Higgs boson hsM can be enhanced by the soft parameter A K and, in order to 
meet the stringent LEP constraints, the hsM tends to decay into the singlet Higgs pairs 
aa or hh instead of 66, which will give a multi-muon signal for hgM produced at the LHC, 
hsM —> aa —> 4/i or h SM — > hh — > 4a — > 8/x. 
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